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Abstract 



A flat membrane with given shape is displayed; two points in the membrane 
are randomly selected; the probability that the separation between the points have 
a specific value is sought. A simple method to evaluate the probability density is 
developed, and is easily extended to spaces with more dimensions. 

1 Introduction 

A most common geometrical problem encountered in exact and natural sciences (engi- 
neering, physics, chemistry, biology, etc) is: given a surface with finite area and definite 
form, and randomly choosing two points in it, find the probability density V{1) that these 
points have a prescribed separation /. 

To see the relevance of the subject consider the following example taken from biology: 
a live membrane is infected at some spots, and the progress of the infection is suspected 
to depend on the mutual separation between the infected points. A knowledge of the 
mean separation Imean, the mean inverse separation {l~^)mean, and the mean squared 
inverse separation {l~^)mean between points in the membrane, all depend on the probability 
density, and are crucial to an analysis of the process. 

In this note we find the functions V{1) for three finite surfaces widely dealt with: the 
circle, the square, and the rectangle. From these examples the method for obtaining 
V{1) for other surfaces is trivially inferred. The method can also be easily extended to 
three-dimensional spaces, such as spherical drops, or parallelepipeds. 



2 Circular membranes 

In a plane disk with diameter 6 two points A, B are randomly chosen. We want the 
probability V{1) dl that the separation between the points lies between / and I + dl. The 
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probability density V{1) has to satisfy the normahzation condition 



V{l)dl 



Our first concern is: having chosen a point A in the {infinite) euchdian plane, next 
randomly choosing another point B also in the plane, we seek the probability that their 
separation have a value /. That probability clearly is proportional to the measure of the 
locus oi B - the circle with center A and radius / - so the locus has measure 2tiI. However, 
since neither the finiteness of the disk nor its shape were taken into account, the analysis 
of the problem is still incomplete. 

We have to examine all possible line segments entirely embedded in the disk. As a 
matter of fact, the symmetry of the disk permits restrict the study to segments aligned 
in just one direction; we choose the vertical direction, for definiteness. 



Figure 1 




In figure 1 we note that the upper tip 5 of a vertical segment with length / has to lie in 
the lens shaped region enclosed by the circles + y'^ = (5/2)^ and + {y — l)^ = (5/2)^. 
The area of this lens is 

1 



S{1) 



(2) 



Writing the probability density as V{1) = k IS{1) and imposing the normalization condi- 
tion (P we find k = 32 /{n6^), and finally 



16/ 



V{1) = — 6\os-\l/6) - iVW^ 

TTO L 



(3) 



A graph of V{1) is given in the figure 2; also a normalized histogram obtained via computer 
simulation is reproduced, to give confidence in the calculations. 



8 P(l) 



Figure 2 
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3 Square membranes 



Similarly as before, we randomly choose two points in a square with side a, and want the 
probability density that their separation be /; the normalization condition now reads 



/ V{l)dl = 1. 
Jo 



(4) 



Again an overall multiplicative factor / is expected in the expression of V{1) , and we are 
next interested in the line segments lying entirely inside the square. 

We initially consider the segments with length I < a; in this case the symmetries of 
the square allow reduce our study to the segments with slope lying between (j) = and 
= 7r/4, as is seen in the first figure 3. 



il<a) 



a - /cosa 



{I >a) 



a - /cosa 



Figure 3 
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^ cos"'fl7/ 

The locus of the upper tip 5 of a segment with slope and length / is a rectangle with 
area 

S{1, (j)) = {a — I cos (j)){a — I sin 0). 
Taking into account all slopes in the range [0, 7r/4] we write 



(5) 



v(i<a) = ki r sa, 

Jo 



(6) 



For line segments with length / > a the minimum slope is cos ^ a/l, as is clear in the 
second figure 3, so now 



V{l>a)=kl S{1, 

J cos ^ a/l 



(7) 



The normalization constant is obtained from 



r V{1 < a)dl + r V{1> a)dl = 1, 

Jo Ja 



and has value k = 8/a'^. We then have V{1) given by the two expressions 

21 r 

V{l<a) = — l'^ -Aal + Txaf 



(8) 



(9) 
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21 r r 
V{l>a) = — AaVl 



4a^ sin ^ a/l - f - na^ - 2a^ 



(10) 



A graph oiV{l) is given in the figure 4; also a normahzed histogram obtained via computer 
simulation is superimposed for comparison. 



aP(l) 



Figure 4 




4 Rectangular membranes 



We assume a rectangle with sides a and 6 < a; to investigate separations / between points 
in the rectangle we need now distinguish three different possibilities, depending on the 
value of / relative to a and h. See figure 5. 



{l<b<a) a - IcosS 



Figure 5 




\b - Isinif) 
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il>a) 
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b 





b - /sin0 



a - /cos0 

\b - /sin0 



sin ^bll (j) sm^l)/l (f) cos'^a/l 

In all cases the angular probability density V{1, 0) is similar as before, 
V{1, (p) = kl{a — I cos (j)){b — I sin 0). 



The symmetries of the rectangle permit restrict the study to segments with slope from 
= to = 7r/2. When I < b < a we integrate (|Tl|) from = to = 7r/2 and 
obtain V{1 < h). When b<l <a the maximum slope is reduced to = sin^^ b/l, and the 
integration gives V{b<l<a). Finally, when l>a>b the slope ranges from = cos~^ a/ 1 
to = sin~^ b/l, and the integration gives V{l>a). 
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The normalization constant k is still unassigned; to fix it we impose the normalization 
condition 



/ V{l<h)dl+ / V{h<l<a)dl+ / V{l>a)dl = l 

Jo Jb Ja 



(12) 



and find k = A/{ab)^] the probability density V{1) is then expressed in the three stages 



Vil<b) 



4/ 



1^2 - (a + b)l + 7rab/2 



(13) 



V{b<l<a) 



Al 



ab sin"^ b/l-al + aVP - b^ - b'^ /2 



(14) 



V{l>a) 
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[a6(sin ^a/l + sin ^b/l)- 



iVP^ + bVP - - / V2 - Trab/2 -{a^ + b^)/2]. 



(15) 



A graph of V{1) is given in the figure 6, drawn for b = a/2; also a normalized histogram 
obtained via computer simulation is given as illustration. 
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Figure 6 
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5 Discussion 



The problem of investigating separations between points in a given space [|l| recently 
received a strong and unexpected stimulus coming from cosmology 0; more specifically, 
from the so-called cosmic crystallography, which aims to unveil the shape of the universe 



19 



The algorithm formerly used in the present note to find the various probability den- 
sities V{1) closely followed that of cosmologists. However, it soon became evident that 
a new approach was imperative to circumvent the long calculations arising from that 
algorithm, when applied to spaces different from balls. 
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The idea of using normahzed probabihty densities V{1) greatly simphfies the task of 
obtaining mean quantities concerning separations; e.g., we have 



Imean = /'^^ / V{1) dl, {l-')mean = l'' V{1) dl, (16) 

JO Jo 
these are particular instances of the general rule 

[f{l)]mean= f{l)V{l)dl. (17) 



In figures 2, 4, and 6 both abscissa and ordinate were chosen dimensionless; by simple 
inspection we then confirm that the corresponding funtions V{1) are indeed normalized. 

In figure 2 we note that the most probable separation in the disk with diameter 5 is 
Imp ~ 0.42(5, corresponding to the value of / where the function 'P{1) of (^) is maximum; 
while the mean separation is slightly greater, Imean ~ 0.455. We also note that the function 



changes curvature at 1/5 = y2/3 ~ 0.82. 

In figure 4 we see that the most probable separation between points in a square 
with side a is Imp ~ 0.48a, while the mean separation is Imean ~ 0.52a. Although both 
the function V{1) in (|^) and (|iy) and its first derivative are continuous at / = a, the 



second derivative is not: in fact, a'^d'^V{l) /dP abruply changes from the finite negative 
value —4 when I = a — e/2 to the diverging positive value ^/\/e when Z = a + e/2. Of 
course the curvature of 'P{1) changes sign at Z = a. We still note in figure 4 the rapidly 
decreasing density of separations when l/a approaches v^; this was already expected, 
since these large separations correspond to segments with both endpoints in diagonally 
opposite corners of the square, and corner regions are small in comparison with the whole 
square. 

In figure 6, corresponding to the rectangle, we again note that both V{1) of (p!^)-(p!5|) 
and its first derivative are continuous throughout < / < ^/a? + W. And again the second 
derivative shows infinite discontinuity, now aX I = h and also at / = a; nevertheless V{1) 
changes curvature only at l = h. When h diminishes relative to a we find that the graph 
of 'P(Z) gradually resembles a right triangle; when b/a ^ the plot is a straight line 
going from (0, 2) to (1, 0), as in the figure 6 of ref. [jl| or figure 3 of ref. |l^: namely, 
aP{l) =2(1 -//a). 

In extending the present note to three-dimensional euclidean spaces one should replace 
the overall multiplicative factor / in the probability density with a factor /^; this is because 
the locus of the points that are at a distance / from a fixed point in three-space is a two- 
dimensional sphere, whose Anl increases with / . 
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